Report  No.  .fr.fr.  .(R£Y-- ). 


t-y 


UNITED  STATES 
DEPARTMENT  OF  THE  INTERIOR 

BUREAU  OF  MINES 
HELIUM  ACTIVITY 

HELIUM  RESEARCH  CENTER 
INTERNAL  REPORT 


DERIVATION  OF  FORMULAS  FOR  EVALUATING  THE  STANDARD  ERRORS  IN  B, 

C,  AND  N  WHICH  APPEAR  IN  THE  FUNDAMENTAL  EQUATION  FOR  THE 

BURNETT  EXPERIMENT:  Z  =  (Z  /PjNr  P 

r  o  o  r 


B.  J.  Dalton 


Robert  E.  Barieau 


BRANCH 

PROJECT 


Branch  of  Fundamental  Research 


no.  _ 


DATE 


July  1965 


AMARILLO,  TEXAS 


' 


' 


*  W?S(dM5 


ib  22009  z&T- 


Report  No.  66  (Rev.) 


HD 

QleUO 

-H43 

M5t 

no. 

66 


HELIUM  RESEARCH  CENTER 
INTERNAL  REPORT 


DERIVATION  OF  FORMULAS  FOR  EVALUATING  THE  STANDARD  ERRORS  IN  B, 
C,  AND  N  WHICH  APPEAR  IN  THE  FUNDAMENTAL  EQUATION  FOR  THE 

BURNETT  EXPERIMENT:  Z  =  (Z  /P  )Nr  P 

r  o  o  r 


By 


B.  J.  Dalton  and  Robert  E.  Barieau 


Branch  of  Fundamental  Research 


Project  4335 

July  1965  uh  UND  MANAGEMENT  LIBRARY 

J  JLDG.  50, 

DENVER  FEDERAL  CENTER 
P.O.  BOX  25047 
DENVER,  COLORADO  80225 


2 


CONTENTS 


Page 

Abstract . . . .  3 

Introduction . .  .  4 


Method  of  obtaining  the  least,  squares  solution 
for  the  constants  appearing  in  the  equation: 

Z  =  (Z  /P  )Nr  P  . . . 

r  o  o  r 

Expressions  for  determining  variances  and 

covariances  of  the  constants  evaluated  .  .  . 

Evaluation  of  the  variance  of  the  P  ,  ,  x ' s 

r (ca 1c) 

and  any  other  calculated  P  that  reduces  F  to  zero.  ....  20 


Expression  for  evaluating  the  variance  of  the 

compressibility  factor  .  ...........  24 

References  .............  .  ......  28 


. 


3 


DERIVATION  OF  FORMULAS  FOR  EVALUATING  THE  STANDARD  ERRORS  IN  B, 
C,  AND  N  WHICH  APPEAR  IN  THE  FUNDAMENTAL  EQUATION  FOR  THE 

BURNETT  EXPERIMENT:  Z  =  (Z  /P  )NrP 

r  o  o  r 


by 


B.  J. 


Dal ton— ^ 


and  Robert  E.  Barieau 


2/ 


ABSTRACT 

The  method  of  treating  a  set  of  isothermally  measured  pressures 
P  ,  P^,  P^ ,  ...  ,  P  ,  for  a  Burnett  experiment,  consists  of  expressing 

the  compressibility  factor  isotherm  of  the  gas  in  terms  of  a  function 
of  either  P  or  p  and  evaluating  the  volume  ratio,  N,  and  the  so-called 
virial  coefficients  by  least  squares  solution.  This  report  gives  the 
method  of  least  squares  and  its  application  to  the  fundamental  equa¬ 
tion  for  the  Burnett  experiment. 

Formulas  are  given  for  evaluating  the  constants  which  appear  in 

the  equation  Z  =(Z/P)NP  together  with  formulas  for  calculating 
n  r  o  o  r 

the  standard  deviation,  of  a  single  Pr(0bs)  >  the  standard  deviation  of 
each  of  the  constants  evaluated,  and  the  standard  deviation  of  the 


compressibility  factor 
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INTRODUCTION 

The  Helium  Research  Center  is  presently  engaged  in  a  critical 
examination  of  all  of  the  thermodynamic  data  on  helium  that,  appear 
in  the  literature.  The  long-range  objective  is  the  development  of 
a  single  equation  of  state  for  helium  that  will  allow  all  of  the 
thermodynamic  properties  to  be  calculated  and  will  reproduce  all  of 
the  data  within  the  accuracy  with  which  the  data  are  known,  A  pre¬ 
liminary  objective  is  to  obtain  the  best  values  of  second  virial 
coefficients  and  to  develop  an  equation  that  will  reproduce  the  best 
values  as  a  function  of  temperature. 

The  method  of  treating  a  set  of  isothermally  measured  pressures, 
Pq ,  P^,  P^ ,  .  ..,  Pr j  for  a  Burnett  experiment,  consists  of  expressing 
Z  in  terms  of  a  function  of  either  P  or  p  and  evaluating  the  con¬ 
stants  by  least  squares  solution,  Associated  with  this  function, 
the  following  errors  are  of  interest:  the  standard  deviation  of  a 

single  P  .  .  .  ;  the.  standard  deviation  of  each  of  the  constants  eval- 
r (obs) 

uated;  and  the  standard  deviation  of  the  function. 

The  purpose,  of  this  report  is  to  present  the  method  of  least 
squares  as  applied  to  the  fundamental  equation  for  the.  Burnett,  experi¬ 


ment  together  with  formulas  which  were  derived  for  evaluating  the 
best  values  for  the  constants  and  the  above-mentioned  errors. 
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METHOD  OF  OBTAINING  THE  LEAST  SQUARES  SOLUTION  FOR  THE  CONSTANTS 
APPEARING  IN  THE  EQUATION:  =  (ZQ/Po)Nr  Pr 

The  fundamental  equation  for  the  Burnett  experiment  is  of  the 

form 

Z  =  (Z  /P  )Nr  P  (1) 

r  o  o  r 

where  N  is  the  volume  ratio  of  the  Burnett  apparatus;  Pq  is  the 

initial  starting  pressure;  r  refers  to  the  expansion  number;  P^  is 

t  h 

the  pressure  after  the  r— —  expansion;  Z^  and  Z^  are  compressibility 

factors  at  P  and  P  ,  respectively. 

r  o  r  J 

There  are  two  series  expansions  which  can  be  used  for  repre¬ 
senting  Z’s:  the  Leiden  expansion  in  powers  of  reciprocal  volume, 
and  the  Berlin  expansion  in  powers  of  pressure.  We  represented  the 
compressibility  factor  isotherm  as 

Z  =  1  +  BP  +  CP2  +  . . .  (2) 

r  r  r 

and 

Z  =  1  +  BP  +  CP2  +  . . .  (3) 

o  o  o 

and  it  was  assumed  that  equations  (2)  and  (3)  could  be  truncated 
after  the  third  virial  coefficient.  The  Berlin  expansion  was  chosen 
to  rep?resent  the  compressibility  factor  isotherm  because  all  of  the 
parameters  appearing  in  equation  (1)  could  be  expressed  in  terms  of 
the  original  observations. 

Let  the  functional  relationship  between  the  variables,  P  and 
r,  involving  the  three  parameters  N,  B,  and  C,  be 
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F  =  F(r,  Pr,  N,  B,C)  =0  (4) 

Now  because  of  random  errors  in  the  observed  pressures,  when  P  ,  ,  . 

r(obs) 

is  substituted  in  the  above  expression,  F  will  not  be  exactly  zero. 

Let  F  be  the  value  of  F  when  the  observed  values  of  r  and  P  are 
r  r 

substituted  in  equation  (4).  Thus, 


F  =  F(r,  P  /  .  ,  N,  B,  C) 

r  %  ’  r (obs)  5  5  ' 


(5) 


Now  we  assume  that  r  is  accurately  known.  Therefore,  equation  (4) 

may  be  solved  for  P  .  .so  that  equation  (4)  is  exactly  satisfied 

r  ^ C  3  -LC / 

Thus  , 


F  =  F(r,  P  x  ,  N,  B ,  C)  =  0 

r (calc) 


(6) 


and  P  .  .  is  the  correct  solution  of  equation  (6) . 

r  JLC  ) 

Now  AP  ,  the  residual  of  P  ,  is  the  difference  between  P  ,  ,  N 
r’  r  r (obs) 

and  P  ,  ,  x  .  This  is  not  the  true  random  error  in  the  observed  P 

r(calc)  r 

because  we  do  not  know  the  true  value  of  P  .  However,  we  can  maxi- 

r 

mize  the  probability  that  AP^_  is  equal  to  the  true  random  error  and 
this  is  just  what  the  principle  of  least  squares  does.  The  princi¬ 
ple  of  least  squares  says  that  we  maximize  the  probability  that  the 


AP^_  ’  s  represent  the  true  random  errors  by  minimizing  the  sum  of  the 


squares  of  the  weighted  residuals.  Thus,  we  should  minimize  the 


* 
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function 


R 


r 

=  '  w  (AP  ) 

Jl  /  4  \  *<■ 

r=1  r(obs) 


-Z 


W 


r=1  Pr(obs)  1 


(7) 


where  AP  is 
r 


AP  -  Y  =  P  (  ,  v  -  P  t  N 
r  i  r(obs)  r(calc) 


(8) 


and  evaluate  B,  C,  and  N  so  that 


(f§)  "  2  J  W 

r>  Pr(obs)>  N>  C  r“l  *r(obs) 


p  ,  .  „ ‘ivaB 


0  (9) 


m  ,  2 

r,^r(obs)»®j^  5 


I\(ob;ii 

r(obs) 


0  (10) 


( 


— )  =  2  >  W  Y  1  — - 

P  R  r  6  prf ohm)  A**  ' 

r  ’  Pr(obs)  *  B’  C  r*l  ^  b  ^ 


0  (11) 


In  equations  (7) ,  (9) ,  (10)  ,  and  (11)  ,  W  is  the  weight 

r (obs) 

to  be  assigned  to  the  observed  P  .  If  the  P  1  s  all  have  the  same 

precision  ipdex,  then  they  will  have  the  same  weight  and  W  =  1 . 

r (obs) 

If  the  Pr’ s  do  not  all  have  the  same  precision  index,  then 


W 


P  .  2 

r(obs)  Sp 


r (obs) 


(12) 


' 
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where  L  is  a  constant  and  S„ 


is  the  variance  of  P 


r (obs) 


W. 


r (obs) 

In  a  particular  problem,  it  may  be  necessary  to  assume 

=  I  in  the.  beginning.  However  ,  if  this  is  done  ,  the 


r (obs) 

residuals,  Y .=  [ P  ,  lS  -P  ,  ,  .1 ,  should  be  examined  to  see  if 

i  r (obs)  r(calc)J’ 

there  is  any  statistical  evidence  for  the  residuals  squared  being  a 
function  of  P^  (0bs)  ”  Any  assumPt;ton  as  to  the  variance  being  a 
function  of  P^  .  ^  ;  can  always  be  checked  by  examining  the  residuals 


In  any  event,  W 


is  not  a  function  of  the  constants  to  be  eval 


r (obs) 


ua  ted 


In  order  to  evaluate  the  constants,  we  need  to  linearize  Y 


l 


with  respect  to  the  undetermined  constants.  A  truncated  Taylor's 
series  expansion  was  used  to  do  this. 

In  a  previous  report  (JL)~^ ,  we.  show  that:  the  linearized  normal 


3/  Underlined  numbers  in  parentheses  refer  to  items  in  the  list  of 
references  at  the  end  of  this  report. 


equations  are.  expressible  as 

a]L  AB  +  b].  AC  +  c][  AN  =  (13) 

t 

a2  AB  +  b2  AC  +  c2  AN  =  m2  ( 3-4) 

a  AB  +  b  AC  +  c  AN  =  m 


(15) 


■ 


. 


af  KJ 
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Equations  (13),  (14),  and  (15)  result  from  expanding  Y^,  (SY  /BB)  , 
(BY  /Sc) ,  and  (BY^/BN)  about  the  approximate  solution  Y° ,  ignoring 
second  and  higher  order  derivatives .  The  linearized  coefficients 
AB ,  AC,  AN  are  defined  as 


AB  =  B  -  B°  (16) 

AC  =  C  ~  C°  (17) 

AN  =  N  -  N°  (18) 

where  B,  G,  and  N  are  our  undetermined  constants,  and  B°,  C°,  N° 
are  approximate  values  for  these  quantities. 

The  a's,  b's,  c's,  and  m's  appearing  in  the  normal  equations 
are  given  to  be  (_1)  . 


r/3Y.»o2 

-  >  W  :j  :  i 

u  p  ,  ,  nLABb  J 

r(obs) 


,B2Y.s  o, 

O  /  L 
*  ’*  ^  ■ 


(19) 


a  „  =  b . 


r 

‘P 

Li  P 


j B Y  .  o  BY  ,so 


r=l 


r (obs) 


l\BB 


J/ 


•  B  Y,.vO*i 

iac“i  '  Y:i  IdBSC/1  j 


(20) 


r 

r  tT 
/  p 


2, 


U 


r  =1 


r (obs) 


BY,  -  o  BY,  o 

1  ;  /  L  V.  O  ;  i  \  a 

+  y  ] 


BB  BN 


8  Y  o 

l  \BBBnJ  j 


(21) 


10 


r 

rrt 

;  w 

i  p 


r=.l 


dY .  o2  c)2Y.  o 

r(obs)V*aC  J  1  Sc2  }  J 


(22) 


and 


r 


■  fdY.  o  BY .  o  dV  o 

b3  ■  c2  -  2  WP ,  h  Aksr)  Ksfj  +  Y?  KFaS/  '  (23^ 

r=1  r(obs) 


r.  r,dY.,o2  S2Y.  o 

2wp  Ibri  +  Y°  (ftv  j 

^  r(obs)  oN 


(24) 


m. 


r 

}  w 


Y.-,  o 
,o  l 


r=l 


P  ,  .  .  i  VdB 
r (obs) 


mr 


r 

rr 


,dY.  o 
o  l  \ 


W  ^  „ 

>.  P  x  ,  N  1.  dC 
r (obs) 

r=l 


(26) 


m. 


r 


•  /  W  x 

L  P  e  ,  N  1  VdN  J> 
r=1  r (obs) 


/dY.*  o 
,o  ./  i  \ 


(27) 


^Yi\° 


,dB 


b  !! 


CdF/dB) 


o 


(dF/dP  ,  .  j 

r(calc) 


o 


(28) 


dY.\  o 
i  ] 

dC  » f 


OF/ao 


o 


(3F/aPr(calc)>° 


(29) 


. 


. 
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BY. ,  o 

_ 1 

BN 


(BF/BN) 


o 


<3F/3Pr(calc)> 


o 


(30) 


,92Y.xO 

(-t) 

SB  ' 


(cS2F/aB2)° 

^F/3Pr(calc))° 


2(5  F/^Pr(calc))0(5F/5B)0 


[(SF/3P. 


.,2 


r  (calc) 


o 


+  QF/as)0  (a2F/apr2(caia)° 


L<3F/3Pr(calc)> 


o 


(31) 


a2Yi,  o 


BC 


2 


o2f/sc2)° 

(3F/SPr(calc)> 


2(a2F/acapr(ealc.))°(aF/ac)0 


(3F/3Pr(calc)> 


+ 


(aF/ac)°  (aW(calc/ 


| (BF/BP 


r (calc) 


*"5  3 

>] 


(32) 


2  2 

B  Y.vO  ,B  Y .  x  o 

_ =  f _ i\  = 

bnbb J  \BBSNy 


(b2f/bnbb)° 


(a  F/aBapr(calc))°(aF/aN)° 


(BF/BP 


r (calc) 


o 


[(3F/3Pr(calc)> 


o 


(BF/BB) °(B2F/BPr (ca lc) BN) °  (BF/BB) °(BF/BN) °(B2F/BP2 (ca lc) )° 


[(3F/3Pr(calc)>°j 


(3F/3Pr(calc)> 


o  I 


(33) 
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,32Y..o  ,32YiX  o 
oNoC/  =  V3cW 


(32F/3N3C)°  <a  F/3C3Pr(caic))0(3F/3N)0 


(SF/3Pr(calc)> 


o 


(SF/SP  (  .  .) 

r(caic) 


o 


(3F/3C)0(32F/3Pr(ealc)3N)°  (3F/3C)°(3F/3N)°(32F/3P2fcale)) 


o 


-,2 


(3F/3pr(caic))  j 


F (3F/3P 


o 


r(calc)  _ 


(34) 


,o 


[P  f  ,  v  -  P  f  ,  v] 
r(obs)  r(calc) 


o 


(35) 


(S1)- 


(32f/3n2->° 

(3F/SPr(calc))0 


2(3ZF/3N3Pr(eale))°(3F/3N)0 


(3F/3P 


r  (calc) 


>°T 


+ 


(3F/3N) 0  (a2F/3p2(cale)) 


,<3F/3Pr(calc)>  J 


(36) 


/a 2yi\°  /^2y±)° 

IdCSB J  =  Vaidci 


(52F/dCSB)°  _  (d  F/5BdPr(calc))  (dF/5C) 


(SF/BPr(calc)> 


o 


(BF/5P  (  v) 
L  r(calc) 


o 


(3F/3B)°(32F/3Pr(ealc)3C)°  (3F/3B)0(3F/3C)°(32F/3P2(ealc:))0 


r 


u(3F/3Pr(calc>) 


o 


--,2 


[(aF/3Pr(caic)>°'. 


od 


■ 
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Now  in  order  to  evaluate  the  solutions  of  our  linearized  normal 
equations,  we  need  values  of  first  and  second  derivatives  of  the 
function,  F, 


F  =  F(r ,  P  ,  .  v  ,  N,  B ,  C)  =  0 

r(calc) 


(6) 


or 


F  = 


Z  ,  .  N  -  (Z  /P  )Nr  P  .  ,  \ 

r(calc)  o  o  r(calc) 


=  0 


f  1  +  BP  +  CP 

=  1  +  BP  ,  .  ,  +  CP  ,  ,  v  1  ° 

r(calc)  r(calc) 


-(■ 


o 


-y 


calc) 


=  0 


These  derivatives  are: 


/'M'l 

Vsb; 


“  Pr (calc)  (1  -  N  > 


/$£) 

sac  / 


=  J  p  ,  ,s  -NP 

r(calc)L  r(calc) 


o] 


=  -  r(Z  /P  )Nr 

o  o 


1 


P 


r (ca lc) 


\ 


dF 


5P 


r (calc) 


) 


=  B  +  2CP  .  ,  N  -  (Z  /P  )N 

r(calc)  o  o 


(38) 


i  .  J  t  0/  •  '  ■ 

i  * 


. 


P,  :  U  .  '  ..  r* 


. 


. 


. 


s2f 


dN 


2/ 


-r(r-l)(Zo/Po)Nr"2  Pr(calc) 


(—r 1  )  ■  2C 

r  (calc) 


CiV 


\dCdBJ 


=  0 


®)  "  -  r  Nr_1  Pr(calc) 


(if)  "  P 


o  r(calc) 


cTF 

VdP  ,  .  NdB/ 

r (ca 1c) 


=  1  -  N 


d  F 


dP  ,  ,  N  dC/ 

r (calc) 


=  2P  ,  .  .  -  Nr  P 

r(calc)  o 


o  F 


dP  ,  .  .dN, 

r (ca 1c) 


=  -  r  (Z  /P  ) N 

o  o 


r-1 


f— a-E - 

\3N9P  .  .  . 

r (calc) 


-  ( 


d  F 


IdP  ,  .  .dN 

r  (calc) 


) 


s2f 


dBdP 


r (calc) 


■ 


/  d  F  "] 

VdP  ,  ,  .  SB,' 


r (calc) 


:  _ B2F 

VdCdP 


r  (calc) 


azF  \ 

dP  /  ,  .dC 

r (ca 1c) 


r  _  .... 
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Now  in  evaluating  the  best  values  for  B,  C,  and  N,  we  must  solve 
the  linearized  normal  equations  by  an  iterative  procedure.  When  we 
have  the  correct  solutions,  equations  (13),  (14),  and  (15)  will  be 
satisfied  exactly.  The  solutions  to  equations  (13),  (14),  and  (15) 
are  (J.) 


where 


D  AB  = 
o 

Dlml  +  D2m2  +  D3m3 

(39) 

D  AC  = 
o 

D,  m.  +  D_m_  +  D^m„ 

4  1  5  2  o  3 

(40) 

D  AN  = 
o 

D^nn  +  Dgm2  +  D^m^ 

(41) 

V 

b2°3  '  b3C2 

(42) 

ii 

D2  =  b3Cl  "  bl°3 

(43) 

D7  " 

D3  =  blC2  "  Vl 

(44) 

D5 

alC3  ‘  a3Cl 

(45) 

« 

00 

II 

D6  =  a2Ci  -  alc2 

(46) 

D9 

=  aLb2  - 

(47) 

D 

o 

D^a^  +  D 2a2  ■*"  ^3a3 

(48) 

'i!  L  ) 


. 
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EXPRESSIONS  FOR  DETERMINING  VARIANCES  AND  COVARIANCES 
OF  THE  CONSTANTS  EVALUATED 


We  now  proceed  to  evaluate  the  variances  of  the  constants  and 
all  of  the  covariances.  To  do  this,  we  apply  the  law  for  the  "Pro¬ 
pagation  of  Errors"  (_3 ,  4)  .  This  law  states  that  if  we  have  a  func¬ 
tion  or  quantity,  say  Q,  that  is  a  function  of  the  independently 
observed  quantities  y  ,  y^ ,  ...  ,  then  the  variance  of  the  quantity 

Q,  is  given  as 


id  J 

i=l  yi(obs) 


(obs) 


(49) 


2  2 

where  S^  is  the  variance  of  Q  and  S  is  the  variance  of  y.,  .  N 

*  y^obs)  l(obs) 

Extracting  the  square  root  of  the  variance,  we  obtain  a  value  on 
the  same  scale  as  the  original  measurements.  This  value,  S^ ,  is 
called  the  standard  error  or  the  standard  deviation  of  Q. 

Now  suppose  we  evaluate  the  variance  of  the  constant  B.  The 
value  of  B  which  we  have  evaluated  is  a  function  of  all  of  the  ob¬ 
served  P  '  s  and  r's.  Since  we  have  assumed  that  r  is  accurately 
r 

known,  then  the  expression  for  the  variance  of  B  is  given  by  the 


equation 


B 


r 

V 

=  / 

r=l 


dB 


2 


dP 


r (obs) 


(50) 


r (obs) 


and  there  will  be  an  equation  similar  to  equation  (50)  for  deter¬ 
mining  the  variance  of  C  and  the  variance  of  N. 

To  evaluate  equation  (50),  we  must  evaluate  (dB/dPr )  for 

,  square  the 


each  j  multiply  this  quantity  by 


r (obs) 


■  I 

.  • 

■  -  : 


- 


' 


. 
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product  and  then  sum  the  product  over  r  from  1  to  r  for  the  r  expan¬ 
sions  . 

In  a  previous  report  (2}  ,  we  have  outlined  the  details  for 
evaluating  the  variances  and  all  of  the  covariances  of  the  constants 
evaluated . 

For  our  particular  problem,  these  var iances- and  covariances  are 
determined  from  the  following  relations  (2): 


B 


D 


o 


2  r 

Di  .1  wp 
r=l 


r (obs) 


2  r  2 

3y.vo  r  Y.vo 

+  d2Zwp  . .  .(A 

r(obs) 
r=l  v  ' 


(— i)' 

\  dB  / 


+  D3iWP 

r=l 


9Y.,o‘ 


r (obs) 


(s1) 


SY.\o,3Y.vo 

i  1  /  i) 


+  2did2  I!  wp  ,  abb  /  \ac 

.  r(obs) 
r=l 


+  2D  D0  ; 

1  3  Zj  P 

r=l 


AY.  ,0/BY.vO 

/  11/  i\ 


r (obs) 


VdB  )  AN  )  +  2D2D3  L  WP 


r=l 


r (obs) 


s'dY  ,\0  ibY .  .  o 

*)  (sr) 


dC 


(51) 


D 


o 


,  A  ,3Y.l02 

D2  \  W  — 1 

4  L  P  ,  .  ,  \dB 
,  r(obs) 
r=l 


+  D5  1  4 


r (obs) 


,3Y..o‘ 

(A) 


,  f  AV° 
+  d6  A  wp  ,  ,  iaiT 

I  u  r(obs) 

r=l 


r=l 

2  r 

+  2D4D5  L  wp  ,  Abb  /  vac  ) 

.  r (obs) 
r=l 


c)Y.xo  AY.,o 

l  \  t  i 


r 

+  2D, )  W 


i  D  >  „ 

4  6  L  P  ,  ,  ABB 

r=1  r (obs) 


BY.* o  BY. n o 

‘Ms1) 


A'1 


BY, ,  o  ,dY ,  o 

+  2d5d6Zwp  lh^)\w) 

,  r(obs) 
r=l 


(52) 
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18 


N 


D 


o 


BY.  o* 
1 


r  2  r 
o  SY.  o  2 

W  — -  +  n  )  W 

7  P  ,  ,  v / SB  7  U8  L  P  ,  ,  ,  SC  . 
r=1  r(obs)  r=1  r(obs) 


r 

e-'-» 


BY  o  BY  o 

1  ll  £  1 


?  r  SY.  o2 

+  d9  l  WP  .  K  +  2D7D8  l  WP  ,  „  aT/  lac 

,  r(obs)  .  r(obs) 

r=l  r=l 


r 

C* 

+  2D  D„  /  W 


BY,  o#?BY.vO 


1D9  L  WP  ,  Abb 

r=1  r(obs) 


i  l 


V 


+  2D  D  W 
•BN  8  9  -  p 

r=l 


r (obs) 


BY.  o  BY.  o 
__i  i  __i  : 

SC  J  \SN  J 


BC 


D 


r  2  r 

BY.,.o  i  SY.  o 

D.D.  )  W„  r— +  D0DC  >  W„  (— ~ 

1  4  ZL  P  /  r  \  SB  /  2  5  L  P  ,  ,  ABC 

r=1  r (obs)  r=1  r  (obs) 


r  >v  2 

C*  ^Yi  ° 

+  D  D  )  W  rji 

3  6£.,  P  ,  ,  N  BN 

..  r(obs) 
r=l 


v-  /SY  o  SY  o 

+  (d2d4  +  DiV  L  WP  ,  „  laF: 

.  r(obs) 
r=l 


r 


+  (d1d6  +  d3d4)  l 


w. 


SY  xo  SY  o 

it)  1st) 


r=l  r(obs) 
r 


.SY.  o  SY  o 

+  <°2D6  +  W  l  WP  ,  ,  (.ar)  la/ 

r=1  r (obs) 


(54) 


(53) 


\ 


2  r  2 

,dY.  o  .  ;SY.vo 

DiD7  L  wp  .  AWJ  +  D2D8  L  WP  ,  AdcT/ 

r=1  r(obs)  r=1  r(obs) 

/V2 

+  D3D9  /,  WP  ,  „ 

,  r(obs) 
r=l 


,3Y  vo  3Y  ,o 

+  (DiDs +  d2d7>  L'h  ,  K  Aar;  \w 

.  r(obs) 
r=l 


r  /SYi\°r3Yi'° 

+  <did9  +  W  Z.wp  ,  ,  Aar)  laiT 

r_^  r(obs) 


^  »3Y..o  BYao 

+  <d2d9  +  D3V  2wp  ,  K  Aar)  (ar) 

r_^  r(obs) 


(55) 


r  2  r  ,  2 

^  Y .  -v  o  ^  /  S  Y  #  ^  o 

V7ZWP  . .  .var)  +  Vs  X  WP  .  „  A ar 

.  r(obs)  .  r(obs) 

r=i  r=l 


BY.s  o' 


+  D  D_  /  W_  n  ^  /i 

6  9  A  P  ✓  ,  N \dN  7 
n  r(obs) 
r=l 


dY.yo,dY^.  o 

/  T  L)  —  Lf  )  /  W_  l  T"~  j  I  — ” 

48  57  /  P  x  ,  N\dB  /  \BC 

,  r(obs) 
r  =  I 


+  (D.  D  +  D  D  )  V  W 


-BY.vO-3Y.-o 

+  (D.  D.  +  D,D_)  )  W  {  — -)  i  —4 

v  4  9  6  iJ  L  p  /  ,  Abb  /  Vsn  / 

r(obs) 


A  y-dY.  oAY.  o 

+  <D5D9  +  D6  V  A  WP  ,  h  ,  \W)  Kl t) 

.  r(obs) 
r=l 


(56) 
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We  have  indicated  that  the  solutions  to  our  linearized  normal 
equations  are  to  be  solved  by  an  iterative  procedure.  When  we  have 
the  correct  solutions,  equations  (13),  (14),  and  (15)  will  be  satis¬ 
fied  exactly.  Therefore,  once  we  have  determined  the  best  values 
for  B,  C,  and  N,  the  remaining  questions  to  be  answered  are:  (1) 
what  is  the  variance  of  the  calculated  P^ ' s  and  any  other  calculated 
P  that  reduces  the  fundamental  equation  for  the  Burnett  experiment 
to  an  exact  identity? ;  and  (2)  what  is  the  variance  of  the  compressi¬ 
bility  factor? 


EVALUATION  OF  THE  VARIANCE  OF  THE  P  ,  .  ' s  AND  ANY  OTHER 

r (calc) 

CALCULATED  P  THAT  REDUCES  F  TO  ZERO 
Because  of  random  errors  in  the  observed  pressures,  equation 
(4)  will  not  be  satisfied  exactly  when  the  observed  P  ' s  are  substi¬ 
tuted  into  this  equation.  Hence,  we  need  to  evaluate  the  variance 

of  the  calculated  P  ' s  which  will  reduce  F  to  zero  for  each  observed 

r 

r  value.  We  do  this  in  the  following  way:  P  .  s  is  a  function 

C  \  C  3.  L  C  ) 

of  the  observed  r's  and,  through  the  constants  evaluated,  is  a 

function  of  all  of  the  observed  P's.  To  evaluate  the  variance  of 

r 

P  .  .  v  ,  we  apply  the  law  for  the  "Propagation  of  Errors"  (_3,  4)  . 

r  C3  ic  ) 

From  equation  (49) ,  we  see  that  this  involves  evaluation  of  the 
quantity 


r (ca lc) 


r  BP  2 

Y  f  r  (ca  lc)  \  g2 

^r(obs)  ^  ^r(obs) 

r=l 


(57) 


Now  in  order  to  evaluate  equation  (57)  ,  we  must  evaluate 

[BP  .  .  n/BP  /  .  ■.  ]  for  each  P  ,  .  ,  ,  multiply  this  quantity  by 
r(calc)  r(obs)  r(obs)  r  n 


-  '  ■  !  A  ; 


»  .  '  r. 
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S  ,  square  the  product  and  then  sum  the  product  over  all  of 

^r (obs) 

the  observed  P  's.  [dP  ,  -  s/dP  ,  .  .]  can  be  determined  from 

r  r(calc)  r(obs) 

equation  (6) .  Suppose  we  differentiate  equation  (6)  with  regard 

to  P  ,  ,  N  ,  holding  r  constant 
r(obs) 


dF 


u 

VdP  f  .  v/  ASP  t  N 

r(caic)  r,B3CPN  r(obs) 


5Fr(calc) '  /dF 
^  +  \dB/ 


SB 


dP  .  .  . 


r ,G pN3Pr (caic)  *r(obs) 


+ 


’2E  :  f - 

*>dC/  ^dP  (  ,  v 

r’B’N’Pr(calc)  r(0bs) 


SC  \  (3FN 

J  VdN ) 


(' 


ML—) 


dP 


r’B’G’Pr(calc)  r(obs) 


0  (58) 


or 


SF  \  / _ r (calc)  |  = 

dP  f  .  J  _  p  AdP  ,  ,  N  / 
r(calc)  r,B3C3N  r(obs) 


(S) 


dB  N 


r  jc  3N 3Pr (calc) 


.dP 


r (obs) 


+ 

\dcJ 


VdN 


dC 

dP  J 

r’B’N’Pr(caic)  r(obs) 


3N 


\  ^  p 

r>BjC’Pr(calc)  r(obs) 


(59) 


and  solving  equation  (59)  for 


[dP  t  .  x/dP  ,  ,  J 
r (calc)  r  (obs) 


we  get 


•  •  ' 

■  -i  ’ 1 


v,  . 


. 
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r (calc) 

r (obs) 


(St) 

V3B/ 


+  1 

Vac/ 


SB 


SP 


r’C>N’Pr(caic)  r(obs) 


SC 


SP 


r’B’N>Pr(calc)  r(obs) 


( 


SN 


SP 


r’BjCjPr(calc)  r(obs) 


SF 


SP 


r (calc) 


r  ,B  3C  3N 


(60) 


Multiplying  equation  (60)  by  S 


r (obs) 


squaring  the  product,  and 


summing  over  all  of  the  observed  P^ ' s ,  we  get 


r (calc) 


,2  / SF' 

’b 


r  ,C  jN  ,P 


+  2S2 

BC  V  SB,' 


/'SF 


i  - 


\sc 


r (calc) 


’C,N’Pr(calc)  r,B’N,Pr(calc) 


,  ,2  2  /3F\ 

+  SC  \SC/  2SBN  SB/ 

r  ,xs  jiN  arr  (calc) 


,  SF 

.SN 


r,C3N,,P  ,  .  r,B3C3P  ,  ,  N 

3  3  3  r (calc)  3  5  r(calc) 


Q2  /MA2 

bN  \SN/ 

NO  r  ,B  3C  ,P 


+  2S2 

+  2bCN  VSC / 


,fSF 


r (calc) 


r ,B ,N 3Pr (calc)  r >B  >G  >Pr (calc) 


SF 


'  2 


SP 


r (calc) 


r  ,B  ,C  3N 


(61) 


■  i 
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from  which  we  can  evaluate  the  (n-i)  values  of  corresponding 

r  (calc) 

to  the  observed  r's. 

Now  we  ask  ourselves  the  question:  how  do  we  calculate  the  var¬ 
iance  of  any  other  calculated  P  that  exactly  satisfies  equation  (6)? 

In  order  to  answer  this  question,  we  must  find  a  value  of  r,  say  r^, 
which  exactly  satisfies  the  equation 


(calc) 


(Z  /P  ) 
v  o  oy 


N 


P 


(calc) 


(62) 


where 


(calc) 


1  +  BP(calc)  +  CP(calc) 


(63) 


Now  suppose  we  have  evaluated  B,  G,  and  N  by  some  means  or  other 
Then  rearranging  equation  (62)  as 


N 


P 


(calc) 


(Z  /P  )  P,  .  . 
o  o  (calc) 


(64) 


and  taking  natural  logarithms  of  equation  (64)  ,  we  get 


In  Z_,  -  In  (Z  /P  )  -  In  P,  .  x 

P ✓  i  \  o  o  (calc) 

_ (calc) 

In  N 


(65) 


from  which  we  can  evaluate  r  ,  and  rp  is  just,  the  expansion  number 

corresponding  to  P,  .  „ 

^  C3 1C J 


.J. 


(1\  nl 
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To  evaluate  the  variance  of  P.  ,  ,  ,  we  must  evaluate  the 

(calc) 


quantity 


(calc) 


t  /3P(calc))2 

L  VdP  .  .  J 

r_jL  r(obs) 


r (obs) 


(66) 


Now  in  order  to  evaluate  equation  (66) ,  we  return  to  equation  (61) 
where  the  terms  involving  derivatives  of  F  are  to  be  evaluated  for 

^(calc)  ’  rP ° 


EXPRESSION  FOR  EVALUATING  THE  VARIANCE  OF  THE  COMPRESSIBILITY  FACTOR 

We  now  proceed  to  evaluate  the  variance  of  Z,  To  do  this  we 

apply  the  law  for  the  "Propagation  of  Errors."  Now  Z,  .  .  is  a 

(  c  a  1  c  ) 

function  of  r^  and,  through  the  constants  evaluated,  is  a  function 
of  all  of  the  original  •^r(0|3S^,s-  From  our  equation, 


'P(calc) 


1  +  BP,  .  +  CP,  . 
(calc)  (calc) 


(63) 


where  is  any  calculated  P  which  exactly  satisfies  equation 

(6) ,  we  see  that  the  variance  in  Z^ 
involves  evaluation  of  the  quantity 


(calc) 


,  as  given  by  equation  (49) , 


SZ 


P  ,  2 

(calc) 


,dP 


(calc) 


r=l 


r (obs) 


r (obs) 


(67) 


. 
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[3Z 


In  order  to  evaluate  equation  (67) ,  we  must  evaluate 

/3P  ,  ,  N  ]  for  each  P  ,  .  .  ,  multiply  this  quantity  by 

P(calc)  r(obs)  r(obs) 

,  square  the  product,  and  then  sum  the  product  over  all  of 


r (obs) 

the  observed  P  '  s,  When  we  do  this,  we  get 

r 


(calc) 


B2  S2  +  4C2P2  .  S2  +  P2  S2 

P(calc)  (caic)  P(calc)  (calc)  B 

+  P^  .  N  S2  +  4BCP  ,  .  N  S2  +  2BP/  ,  N  „ 

(caic)  C  (caic)  P,  .  N  (calc)  B,P,  ,  N 

(calc)  3  (calc) 


2  2  2  2 

+  2BP .  i  \  S_,  _  +  4CP,  .  ,  _ 

(caic)  C5P(calc)  (calc)  B,P(calc) 

+  4CP(calO  .  ,  +  2P(calc)  4 

v,c  a  lc_) 


(68) 


2  2 

The  terms  S  and  S  ,  which  appear  in  equation 

3  (caic)  L^(calc) 

(68)  ,  are  defined  as 


B,P 


(calc) 


3P 

x3P 


r=l 


r(obs)  ^Pr(obs)  Pr(obs) 


(69) 


C,P 


(calc) 


r=l 


Vale)  .  _3C _  „2 

^Pr(obs)  ^Pr(obs)y  Pr(obs) 


(70) 


J  ■■•'6  W  3.' 


II  H  II  n 
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and  these  quantities  are  to  be  evaluated  from  equations  (69)  and 

(70),  which  we  proceed  to  do. 

[BP,  n  x/BP  /  ,  \]  can  be  determined  from  equation  (60), 

(.calc;  r  (.obs; 

where  the  terms  involving  derivatives  of  F  are  to  be  evaluated  for 


2 


P,  i  \  ?  r  •  Multiplying  equation  (60)  by  (3B/BP  ,  ,  .)  S_ 

\C3. 1C )  Jr  IT^ODSyx 

and  summing  over  all  of  the  observed  P^ ! s ,  we  get 


r (obs) 


BP 


(calc)  Y _ oB_ 


r)P  /\B P  P 

r(obs)  r(obs)  r(obs) 


=  S' 


BSP,  .  v 

(caic) 


m  s2  +  ,’M ,  s2 

3B/rD,C,N,P.  ,  .  B  SC  r„,B,N„P,  ,  .  BC 

P  '  (calc)  P’  5  -  (calc) 


iW 


rp,B 2C  5P 


BN 


(calc) 


BF 


BP 


(calc)  rp,B5C3N 


(71) 


Multiplying  equation  (60)  by  (BC/BP  ,  *)  S 

r \ODSj  r 

summing  over  all  of  the  observed  P^'s,  we  get 


and  then 


r  (obs) 


( 


,-BP 

BP 


)  s2 


.  1 


r(obs)'  ^^(obs)  'Pr(obs) 


=  S 


C  jP/  ,  V 
3  (calc) 


d2  +  3E 

^B/  BC  +  V3C 

P’  ’  ’  (calc) 


r  BNP 

rp^^Ycaic) 


4~ 


/B  F 


r  R  r  P 

rp ,D ^ (calc) 


CN 


(72) 


/  BF 

\BP,  ,  J 


(calc)  rp ,B  ,C  3N 


■  ' ' '  1  . 


a.  (.->3 
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Therefore,  the  use  of  equations  (7.1)  and  (72)  in 

2 


will  enable  the  variance  of  a  calculated  Z  ,  S 

Jl  Zj 

mined . 


(calc) 


equation  (68) 
to  be  deter- 


28 
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